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ABSTRACT 

The  prediction  of  dipole  sound  from  the  diffraction  of  turbulence  by  the 
leading  edge  of  a  thick  foil  is  made  with  the  Kirchhoff  integral  for  rigid  surface 
scattering  of  the  stagnation  enthalpy.  The  incident  field  is  determined  from  a 
volume  integral  and  the  rigid  plane  Green  function  using  an  equivalent  form  of 
Howe’s  (1975)  acoustic  analogy  that  is  derived  in  terms  of  the  mean  free  stream 
velocity  and  the  fluctuating  up-wash  velocity  found  in  Sears’(1941)  analysis.  A 
comparison  of  the  measured  and  the  predicted  dipole  sound  made  with  the  foil 
geometry  and  measured  turbulence  statistics  from  Paterson  and  Amiet  (1976) 
shows  good  agreement.  The  thickness  effect  is  incorporated  in  the  governing 
Green  function  for  the  foil  [Howe  (1998a,  2001a)]  and  serves  to  exponentially 
attenuate  the  dipole  sound  pressure  spectrum  by  the  product  of  the  convection 
wave  number  and  half  the  maximum  section  thickness.  The  dipole  sound  from  a 
foil  cutting  through  a  mean  shear  layer  is  then  calculated  using  an  acoustic 
analogy  where  the  source  has  retained  the  mean  shear  term.  The  ratio  of  the 
dipole  sound  of  the  mean  shear  source  to  the  source  without  mean  shear  from  the 
earlier  calculation  is  determined  to  be  proportional  to  the  ratio  of  the  mean  shear 
to  the  frequency  of  the  sources.  Estimates  of  the  dipole  sound  with  and  without 
the  mean  shear  source  are  made  for  Olsen  and  Wagner’s  (1982)  experiment. 
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1.  Introduction 


An  airfoil  in  a  turbulent  flow  generates  a  net  distribution  of  acoustic  dipoles  primarily  at 
the  leading  and  trailing  edges.  In  these  regions,  the  dipole  surface  normal  stresses  induced  by  the 
turbulence  are  not  canceled  by  their  images  due  to  the  surface  curvature  [Powell  (1960), 
Meecham  (1965)].  Sears  (1941)  modeled  the  incompressible  unsteady  lifting  force  associated 
with  the  dipoles  by  representing  the  foil  as  an  acoustically  compact  two-dimensional  strip.  His 
analysis  included  a  trailing  edge  Kutta  condition  and  showed  for  fluctuating  velocities  much 
smaller  than  the  free  stream  velocity  that  the  unsteady  lifting  force  was  proportional  to  the 
product  of  the  mean  firee  stream  velocity  and  the  fluctuating  “up- wash”  velocity.  The  application 
of  the  Kutta  condition  essentially  removed  this  edge’s  contribution  to  the  incompressible  part  of 
the  total  unsteady  lift  for  reduced  frequencies  greater  than  one  (based  on  the  chord  length) 

[Howe  (2001b)]. 

When  the  airfoil  chord  length  is  acoustically  non-compact,  the  dipole  sound  produced 
from  the  diffraction  of  the  turbulence  by  either  the  leading  or  trailing  edge  is  subsequently 
“back-scattered”  by  the  opposite  edge  of  the  foil.  Landhal  (1961),  Adamczyk  (1974),  Amiet 
(1975),  and  Martinez  and  Widnall  (1980)  accounted  for  the  back  scattering  by  the  trailing  edge 
of  the  leading  edge  dipole  sound  for  a  surface  without  thickness.  Roger,  Moreau  and  Wang 
(2002)  extended  Amiet’ s  (1975)  solution  with  a  second  scattering  correction  term  for  application 
to  a  trailing  edge  noise  problem.  Howe  (2001a)  constmcted  a  Green  function  to  account  for  the 
multiple  back  scattering  by  both  edges.  The  Green  function  also  included  the  dependence  of  the 
shape  of  the  edge  for  the  diffraction  of  the  incident  turbulence  [Howe  (1975)].  In  practice,  the 
use  of  a  Green  function  to  predict  the  dipole  sound  produced  by  the  diffraction  of  a  turbulent 
flow  requires  that  the  spatial  distribution  of  the  incident  sources  be  prescribed  and  remain  frozen 
within  the  diffraction  zone.  This  “rapid  distortion”  approximation  attributes  the  diffraction  to  the 
Green  function  and  does  not  account  for  any  interaction  between  the  incident  sources  and  their 
images  [Howe  (1999)].  This  approximation  is  satisfactory  as  long  as  the  turbulence  fluctuations 
are  less  than  10%  of  the  free  stream  velocity  [Grace  (2001)]. 

The  acoustic  analogy  of  Howe  (1975)  has  been  utilized  with  a  shape  dependent  Green 
function  to  estimate  the  dipole  sound  from  the  diffraction  of  the  sources  within  turbulent 
boundary  layers  by  trailing  edges  of  various  contours  [Howe  (1998a  and  b),  (1988),  (2000)]. 
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Calculations  of  the  shape  dependence  of  the  unsteady  lift  for  turbulence  encountering  the  leading 
edge  have  been  restricted  to  numerical  models  with  idealized  representations  of  the  flow  acoustic 
source.  The  spatial  domain  formulation  of  the  blade  vortex  interaction  (B  VI)  models  of 
Martinez  and  Rudzynsky  (1997)  and  Grace  (2001)  predict  a  reduction  in  unsteady  lift  with 
increasing  frequency  for  a  thick  foil  relative  to  a  foil  without  thickness.  These  trends  agree  with 
the  leading  edge  dipole  sound  measurements  made  by  Olsen  and  Wagner  (1982)  with  airfoils  of 
varying  thickness.  The  BVI  calculations  are  also  substantiated  by  the  difference  between  the 
measured  dipole  sound  spectra  from  a  NACA  0012  airfoil  and  the  sound  spectra  predicted  by  a 
zero  thickness  theory  [Paterson  and  Amiet  (1976, 1977)]  as  illustrated  in  Fig.  1.1.  Paterson  and 
Amiet  ascribed  the  discrepancy  between  their  theory  and  data  to  the  thickness  of  the  foil  for 
convection  wavelengths  of  the  incident  turbulence  that  were  smaller  than  the  foil  thickness. 

Grace  (2001)  showed  that  calculations  performed  in  the  wave  number  domain  of  the  unsteady  lift 
response  of  thick  foils  to  “gusts”  incident  to  the  leading  edge  do  not  capture  the  observed 
thickness  dependence. 

The  current  work  uses  an  acoustic  analogy  and  a  shape  dependent  Green  function  to 
predict  the  lifting  dipole  sound  from  the  leading  edge  diffraction  of  homogeneous  turbulence  by 
a  thick,  symmetric  foil  at  a  zero  angle  of  attack  to  the  mean  flow.  The  foil  thickness,  h,  is 
restricted  to  being  acoustically  compact,  QiK«V)  while  the  chord  length,  c,  can  be  greater  than 

an  acoustic  wavelength,  X,=  ,  (cfco>l). 

In  Section  2,  the  portion  of  the  Howe  (1975)  aero-acoustic  source  that  generates  the 
lifting  component  of  the  dipole  sound  is  expressed  in  terms  of  the  mean  free  stream  velocity  and 
the  fluctuating  up-wash  velocity  that  also  arises  in  Sears’  (1941)  analysis.  The  source  is  tailored 
to  be  linear  in  fluctuating  variables,  without  mean  shear  [Blake  (1986),  Howe  (1989)]  and 
spatially  homogeneous.  The  near  field  of  the  quadrupole  radiation  on  a  rigid  plane  is  then 
obtained  fi-om  the  volume  integral  of  this  source  with  the  appropriate  Green  function. 
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SPECTRUM  LEVEL.  dB 


In  Section  3,  the  leading  edge  lifting  dipole  sound  produced  by  a  foil  is  modeled  by  the 
Kirchhoff  integral  for  rigid  surface  diffraction  of  the  stagnation  enthalpy.  The  Green  function 
incorporates  the  effect  of  the  foil  thickness  [Howe  (1998a)]  and  the  acoustic  non-compactness  of 
the  airfoil’s  chord  length  [Howe  (2001a)].  The  incident  surface  stagnation  enthalpy  is  expressed 
in  terms  of  the  mean  velocity  and  the  fluctuating  up-wash  velocity  spectrum  obtained  from 
Section  2.  This  formula,  along  with  the  measurements  of  the  incident  flow  turbulence  reported 
by  Paterson  and  Amiet  (1976, 1977),  are  used  to  predict  the  measured  dipole  sound  for  a 
NACA0012  airfoil  (/i/c=0.12)  and  the  corresponding  zero  thickness  theory  (/i/c=0)  for  a  variety 
of  reported  flow  speeds. 

In  Section  4,  the  acoustic  source  is  tailored  to  be  linear  with  respect  to  fluctuating  terms 
and  to  contain  the  mean  shear  [Kraichnan  (1956),  Chase  (1980),  Blake  (1986)].  This  source  is 
specified  to  be  homogeneous,  and  the  imposed  surface  pressure  spectrum  from  the  quadrupole 
radiation  on  a  rigid  plane  surface  is  obtained.  The  calculated  surface  pressure  spectrum  for  the 
mean  shear  flow,  when  combined  with  the  results  of  Section  2,  yield  a  prediction  formula  for  the 
leading  edge  lifting  dipole  sound  of  a  thick,  acoustically  non-compact  foil  cutting  through  a 
shear  layer.  Estimates  of  the  measured  leading  edge  dipole  sound  by  Olsen  and  Wagner  ( 1982) 
are  made  for  incident  sources  with  and  without  the  mean  shear  term.  These  estimates  establish 
the  frequency  below  which  the  mean  shear  source  dominates  the  source  without  mean  shear. 
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2.  Wall  Stagnation  Enthalpy  Spectrum  from  Homogeneous  Turbulence  Without  Mean 
Shear 

The  stagnation  enthalpy  spectrum,  B,  imposed  on  a  rigid  plane  (y„  =  0,  yj)  will  be 

calculated  from  the  quadnipole  radiation  of  a  spatially  homogeneous  turbulence  field  as 
illustrated  in  Fig.  2.1.  The  turbulent  flow  is  specified  to  have  no  mean  shear  so  that  the  resultant 
wall  stagnation  enthalpy  spectrum  will  serve  as  the  incident  field  that  will  be  used  in  Section  3  to 
calculate  the  leading  edge  noise  from  homogeneous  turbulence  in  the  absence  of  mean  shear. 

The  radiating  stagnation  enthalpy,  B,  is  described  by  the  acoustic  analogy  of  Howe  (1975)  for 
vortex  sound 

(2.1) 

where  and  U  represent  respectively,  the  sum  of  the  mean  and  fluctuating  vorticity  and 
velocity.  The  stagnation  enthalpy  is  expressed  in  terms  of  the  acoustic  pressure,/?,  the  fluid 
density,  p,  and  the  fluctuating  velocity,  m: 

B=  Jrfp/p  +  O.SM^  .  (2.2) 

The  component  of  the  hydro-acoustic  source  that  contributes  to  the  lifting  dipoles  with 
axes  aligned  normal  to  the  free  streamlines  [Blake  (1986),  Howe  (1989)]  shall  be  used  to  model 
the  incident  source 

V-(a’xt/>i-(n'3xC/,)(y,0.  (2.3) 

The  primed  (')  variables  are  in  the  time  domain,  and  the  subscripts  1, 2  and  3  denote 
respectively  the  streamwise,  spanwise,  and  surface  normal  directions.  The  source  is  linear  with 
respect  to  the  fluctuating  terms,  and  only  the  spanwise  vorticity,  Q3,  contains  the  fluctuating 


Homogeneous  turbulence  over  a  rigid  plane  and  imposed 
wall  stagnation  enthalpy,  B(yi,  0,  y^). 


components  of  the  source.  The  subsequent  planar  wave  number  spectrum  of  the  near  field 
radiation  is  given  by  Blake  ( 1 986) 


(2.4) 


where  the  unprimed  variables  have  been  Fourier  transformed  from  time  to  frequency  and 

^1,3  =  (^i>  ^3)*  -^1.3  =  (^i»  .^3)*  ^ Xi,j=dxidxy  When  the  in-plane  Fourier  transform  of  the  Green 
function 


r  rik) 


'*13*13^2^  _;C _ -lKryi.1 


(2.5) 


is  substituted  into  Eq.  (2.4),  the  stagnation  enthalpy  wave  number  spectrum  becomes 

1  fff  <9  Jyii'W 

JJJ/^(>')^(«3  , 

(2;r)  ^2  7(k) 


(2.6) 


where 


and  =colc^ ,  is  the  acoustic  wave  number.  Rearranging  the  terms  of  Eq.  (2.6), 


B{k,a))  = 


(2Ky 


y(k) 


(2.7) 


while  taking  the  Fourier  transform  of  the  source  variables  with  respect  to  the  in-plane 
coordinates  (denoted  by  the  over  bar): 
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■foo 

^  {271) 


and  upon  substitution  of  Eq.  (2.8)  into  Eq.  (2.7),  in  the  absence  of  mean  shear,  yields  the 
following  expression  for  the  wall  stagnation  enthalpy  wave  number  spectrum: 


j.  3 -  g'WJ 

B{k,a))=  J  dy2Ui-^{fl^{kii,y2,ci)fji 


Howe’s  (1975)  acoustic  source  in  Eq.  (2.9)  can  be  expressed  in  terms  of  the  fluctuating 
up-wash  velocity,  Mj,  and  free  stream  velocity,  (/, .  The  solenoidal  velocity  wave  number 
spectrum  can  be  expressed  in  terms  of  the  vorticity  wave  number  spectrum  [Howe  (1998a]: 


u(k)  =  ik  X  Q(k)  /\k\ 

I  k  P=  k^  -hk^  "^k^  • 


(2.10) 


The  up- wash  velocity  wave  number  spectrum  becomes: 


U2(k)=  i(k^xC^J(k))2/\kf  . 


(2.11) 


Rearranging  Eq.  (2. 1 1), 


■U,U2(k)  =  -i{Cl2(k)xU,)2 


(2.12) 


and  taking  the  inverse  Fourier  transform  Eq.  (2.12)  with  respect  to  k2. 


7  I  ifc  P  - 


(2.13) 


9 


and  then  substituting  Eq.  (2.13)  into  Eq.  (2.9),  yields  the  wall  stagnation  enthalpy  wave  number 
spectrum; 


0  »V2  —  *1  yik) 


(2.14) 


In  the  absence  of  mean  shear,  Eq.  (2.14)  becomes; 


t  giyim  f  I L  |2 

Bik,Q))=  J  dy2i——{  J  dk^{ik^y’^^^ {i—U ^U2iK,k^,k^,(o))) 
0  yw 


(2.15) 


The  complex  conjugate  (*)  of  the  wall  stagnation  enthalpy. 


j  I  d):', (-ft', 

0  71*1  *1  U-lo) 


with  Eq.  (2.15)  enables  the  wall  stagnation  enthalpy  wave  number  spectrum  to  be  determined 


(k,Q))S(k  -k')Sio)-0)')=<  Bik,a))B*  (*',©' )  > 


(2.17) 


where  the  angle  parentheses  (<  >)  signify  an  ensemble  average.  Substituting  Eqs.  (2.15)  and 
(2.16)  into  (2.17)  gives  the  wall  stagnation  enthalpy  wave  number  spectmm 


+tM  +00 


J  dy^fdy'^lj  +  k]  +  klfQ>  (k, ,  3,0)]  (2.18) 

0  0^  ^  I  y{k)  I 


where  the  up-wash  power  spectrum  is  also  obtained  as  in  Eq.  (2. 17); 


^u^u^(k,a))5ik-k  )6(co -(!))=<  u^(Jc,(o)u2  {Jc ,oi) >  . 


(2.19) 


In  order  to  facilitate  the  integration  of  Eq.  (2.18),  the  up-wash  wave  number  frequency  spectrum 
is  expressed  in  a  separable  form, 

where  all  of  the  wave  number  spectra  are  normalized  accordingly: 

+0O 

ldkj<l>j(kj)  =  l.  (2.21) 


The  streamwise  (1)  and  vertical  (2)  wave  number  spectra  will  use  Corcos’(1963)  representation: 


m)= 


(2.22) 


Referring  to  Fig.  2.1,  the  convection  wave  number  Xc^=  -co/Ui  for  gi>0.  The  sign  of  the 
convection  velocity  has  no  bearing  on  the  value  of  the  radiated  noise  power  spectrum.  This  will 
become  apparent  since  all  of  the  subsequent  expressions  for  the  incident  surface  pressure 
spectrum  involve  even  powers  of  the  streamwise  wave  number,  fc,.  The  spanwise  (3)  wave 
number  spectrum  will  be  approximated  with  a  Gaussian  form  used  by  Blake  (1971): 


^3  (^3)  ~  h  2^^  ’  ^3  “  1  XO)  . 


(2.23) 


Substituting  Eqs.  (2.22)  into  Eq.  (2.18)  and  integrating  Eq.  (2.18)  with  respect  to  kj  yields: 

Y  *C  .  g‘y2Y(kyiy'2r\k)  j 

~  J  ^3^2  J  ^yi^l  (^,3’®)]  fy(^)P 

2(-l+Xc^l^  +  klJ^fnSiy^  -  y2)X^  -2(-H-  +  k^l^)7cS\y,  -  y,)^ 

+27c5‘^iy,-y,)t;^]} 
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where  5"  is  the  nth  order  derivative  of  the  delta  function.  The  relationship  between  the 
convolution  of  the  nth  order  derivative  of  the  delta  function,  with  a  function,/,  and  the  nth  order 
derivative  of/is  given  in  Eq.  (2.25)  [Lighthill  (1958)]; 


idyJ(y2)S''iy2-y2)  =  i-iyf''(y2)  OKy'^Koc .  (2.25) 

0 

If  the  integral  scale  is  much  smaller  than  the  spatial  extent  of  the  turbulence,  then  the  vertical  (2) 
correlation  function  can  be  represented  as  a  product  of  the  integral  scale  and  the  delta  function: 

g-\y2-yi\ii2  =  45(^2  -  yj) ;  4  « 1 .4t/,  /ti)  =  1 .4*;*  .  (2.26) 

Substituting  Eq.  (2.26)  into  the  expression  for  the  wall  stagnation  enthalpy  wave  number 
spectrum  [Eq.  (2.24)]  yields; 


f  f  ,  ,  1 

J  dy:,]dy^U,^  ik^^,o})]  .  .2  ' 

00  l7(*)r  k{jt 

{l2[t^\5{y2-y2)l2(-l+2k^ll  +  2klll-k^X -k*l*)K)+  (2.27) 

2(-l+  ^'4^  +  kll^y7uS(y2-y'2)t^ -2(-l+ k^ll  +  kl^)itS{y2-y\)ty 

+2K8‘\y^-y\)r^])  . 


The  forms  of  the  wave  number  spectra  in  Eqs.  (2.22)  and  (2.23)  yield  peak  values  for  the 
streamwise  and  spanwise  wave  number  spectra  at  it,=  k^  =-,+oyf/^  and  ^^3=0,  respectively.  For 
k\>k^k-i,  the  exponential  phasor  in  Eq.  (2.27)  may  be  approximated  as: 

m^ik, . 

Integrating  Eq.  (2.27)  with  respect  to  >2. 
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(2.28) 


7  1 

{ hi  li  ih  (-1 + 2fef +  2klll  -  k^h  -  '2-klklh  -  kt  h  )^) + 

2i-l+klll  +  klll)^M^-2i-l+k^ll+klll)7c4k^l^+^^^  . 


and  then  integrating  Eq.  (2.28)  with  respect  to  yields  the  wall  stagnation  enthalpy  frequency 
wave  number  spectrum: 


{ hi  li  ih  (-1  +  2k^ll  +  ^kHl  -  kth  - 2k\klh  -  ktlt  )^)  +  (2-29) 

2(-l+kfll  +klllf7d^-2(rl+kfll  +  k^ll)7c4k^l^  +2Kl6k*lf]}  . 

The  wall  stagnation  enthalpy  frequency  spectrum  is  determined  by  integrating  Eq.  (2.29) 
with  respect  to  ki  and  k^.  Using  approximate  forms  of  the  wave  number  spectra  of  the  up-wash 
velocity  facilitates  the  integration.  Accordingly,  since  the  streamwise  wave  number  spectrum  of 
the  up-wash  velocity  [Eq.  (2.22)]  is  peaked  at  the  convection  ridge,  and  since  the  rest  of  the 
integrand,  denoted  by/Ck^ks),  of  Eq.  (2.29)  is  varying  slowly  in  this  region,  the  streamwise  wave 
number  spectrum  shall  be  represented  as  a  delta  function.  Thus: 


+oa  +00 

i dk,f(k,.k,)<h(ic,)‘  I ‘iKfiK-ih)S{i‘,-K)= f{KA)  ■ 


(2.30) 


The  wall  stagnation  enthalpy  becomes: 


/g-(t3'3/2)^  1  I  1 

.«>)  =  £', ’  2^  kl  +  kl^'W, 

{/,K’(4(-1  +  le^li  +  2kll}-kX  - 2kll4li  - kX)it)  + 

2(-l+ -2(-l+ +  +  ■ 


(2.31) 
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Integration  of  the  wall  stagnation  enthalpy  wave  number  spectrum  with  respect  to  yields  the 
wall  point  stagnation  enthalpy  frequency  spectrum.  This  involves  three  types  of  integrals, 
/i=0,2,4,  where 


(2.32) 


Approximating  the  error  function,  Erfix)~x  for  x<l,  gives  the  first  integral,  4  [n=0;  Eq.  (2.32)] 


/o  =  Erf[l/2-^ijk^  12  =  0.6 


(2.33) 


The  second  integral,  /j  [n=2;  Eq.  (2.32)]  yields: 


'1  = 


(2.34) 


Numerically  evaluating  /j  gives: 


,  _  2-1. 12-1.36(1. 12)' V^+1.36(1.12)'V^  0.56  . 
h - - OAlc,. 


(2.35) 


The  third  integral,  [n=4;  Eq.  (2.32)],  becomes: 


-/jM + Ik^ll  -  ^^^KEIfO.  )] 


(2.36) 


This  simplifies  to 
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Taking 


^3  =  1.12, 

/  _zl:l^[_4  +  2(1.12)^-1.36(l.l2)^V^+1.3e(L12)^V^-0.56] , 

2(1.12) 

/.  =  1.19A:^ 


(2.38) 


Substituting  the  three  integrals,  /q,  h,  h  into  the  expression  for  the  point  stagnation  enthalpy 


gives 


£/?<>,„(<»)— ((j-^[-0-6+2(1-12)"»«+2(0.4)(1.12)=-(1.12)'o.6 

-2(1. 12)"  0.4-1. 19(1. 12)'‘];r) 

+-^^[-0.6+ 0.6(1.12)^ +  0.4(1.12)^f 


(2.39) 


V  L  I  .i-A-y  ■  w.  /  J 

(1.12)‘ 

.-^5^[-0.6+0.6  {1.12)"+0.4  (1.12)"]+5£l^}  . 
(1.12)  (1-12) 
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The  component  of  the  wall  stagnation  enthalpy  power  spectrum  that  contributes  to  the  hfting 
dipole  sound  expressed  in  terms  of  the  mean  flow  velocity  and  the  fluctuating  up-wash  power 
spectrum  becomes: 


(2.40) 


3.  Dipole  Sound  from  a  Thick,  Symmetric  Foil  in  a  Turbulent  Flow  Without  Mean  Shear. 

In  Section  2,  the  surface  stagnation  enthalpy,  B,  was  calculated  on  a  rigid  plane  due  to  the 
near  field  acoustic  radiation  from  the  component  of  Howe’s  (1975)  acoustic  source  that  gives 
rise  to  the  lifting  components  of  the  dipoles.  The  source  was  specified  to  be  without  mean  shear 
and  was  expressed  in  terms  of  the  firee  stream  velocity  and  the  fluctuating  up-wash  velocity.  In 
this  section,  the  dipole  sound,  .  due  to  the  diffraction  of  these  sources  by  the  leading  edge  of 
a  thick,  symmetric  foil  will  be  calculated  using  a  rigid  surface  Green  function,  G(x,y,CO)  and 
incident  “blocked”  surface  stagnation  enthalpy,  B ,  that  satisfies  the  Kirchhoff  equation 


B,,j(x,co)  = 


j_  ff  dB{y„y^  =0,y3»^) 

2 


Gix,y,ct))dS(y) , 


(3.1) 


where  B,^  is  the  far  field  radiated  stagnation  enthalpy  and  y„  is  normal  to  the  surface  and  is 
directed  into  the  fluid.  At  low  flow  Mach  numbers.  Me  =Ulco<l,  the  radiated  stagnation 
enthalpy,  is  related  to  the  acoustic  pressure  in  the  far  field,  [Howe  (1978)] 

P..=5™,-p/(l+M,(x, //?)),  (3.2) 


where  MXxJR)=Mor  is  the  flow  Mach  number  component  in  the  observer  direction.  The 
acoustic  pressure  can  then  be  determined  in  the  far  field  from  Eqs.  (3.1)  and  (3.2): 


2{\+MeixJR)) 


The  coordinates  for  the  radiation  field  are  illustrated  for  a  foil  of  thickness,  h,  chord 
length,  c,  and  span,  L,  in  Fig.  3.1.  The  foil  encounters  homogeneous  turbulence  with  convection 
velocity,  The  cross  section  of  the  foil  [Fig.  3.2(a)]  shows  that  the  turbulence  flows  past 

both  sides  of  the  leading  edge  of  the  foil.  The  half-plane  “slab”  scattering  geometry  that  is  used 
to  model  the  leading  edge  diffraction  of  the  incident  stagnation  enthalpy  [Eq.  (3.3)]  is  given  in 
Fig.  3.2(b).  The  incident  stagnation  enthalpy  is  imposed  by  the  turbulence  on  both  sides  of  the 
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Figure  3. 1  Coordinates  for  the  acoustic  field  of  the  NACA  0012  airfoil  with 

Span,  L,  thickness,  h,  and  chord  length,  c,  and  convection  velocity,  17^=- 
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leading  edge  for  leading  edge  force  estimate  from  the  diffraction  of  the  incident 
surface  stagnation  enthalpy,  B. 


foil.  A  leading  edge  Kutta  condition  is  not  applied  and  the  vortical  sources  are  prescribed  to 
follow  the  potential  flow  streamlines  associated  with  the  foil  geometry  and  the  free  stream 
velocity  indicated  in  the  figures.  The  surface  of  the  foil  is  defined  to  have  no  discontinuities  in 
the  surface  curvature.  Even  though  the  sources  accelerate  due  to  the  curvature  of  the  potential 
streamlines  in  the  vicinity  of  the  leading  edge,  it  will  still  be  assumed  that  the  incident  source 
wave  number  spectmm  is  peaked  at  the  convection  ridge  associated  with  the  velocity  upstream 
of  the  leading  edge. 

The  dipole  sound  spectrum  due  to  the  diffraction  of  the  incident  turbulence  by  the  leading 
edge  shall  be  modeled  using  a  surface  pressure  scattering  formulation  with  Howe’s  (1998a  and 
b)  trailing  edge  Green  function  for  a  semi-infinite  hard  surface  with  thickness,  h  and  rounded 
edge. 


pp  y  ^  3\  j 


(3.4) 


where  ^pp  (x,co)  and  (<u)  denote  respectively  the  radiated  acoustic  pressure  spectrum  and 
the  incident  surface  pressure  spectrum.  The  right  hand  side  of  the  equation  has  been  multiplied 
by  a  factor  of  two  in  order  to  account  for  the  two-sided  turbulent  flow  around  the  airfoil’s 
leading  edge.  Referring  to  Eq.  (2.27),  incident  sources  at  the  convection  wave  number, 

=  (OfUc,  and  location,  y^,  normal  to  the  surface  have  weighted  contributions  to  the  surface 


pressure  spectrum  that  are  exponentially  attenuated  by  their  product: 


(3.5) 


The  thickness  correction,  in  Eq.  (3.4), 


suggests  that  the  sound  generated  by  a  foil  with  thickness,  h,  in  a  spatially  homogeneous 
turbulence  field  can  be  alternatively  represented  by  a  foil  with  zero  thickness  where  the  incident 
sources  are  off-set  from  the  surface  in  the  surface  normal  direction  by  a  distance  equal  to  one 
quarter  of  the  maximum  section  thickness,  lyjb  h/4. 
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Substituting  Eq.  (2.40)  that  relates  the  mean  velocity  and  upwash  velocity  to  the  wall 
pressure  stagnation  enthalpy  and  Eq.  (3.5)  that  gives  the  Mach  number  correction  into  Eq.  (3.6) 
yields  the  leading  edge  dipole  sound  pressure  spectrum  due  to  the  diffraction  of  the 
homogeneous  turbulence: 


(3.7) 


The  multiple  “back-scattering”  of  the  acoustic  waves  generated  by  the  leading  edge 
diffraction  of  the  turbulence  by  the  edges  of  the  foil  occur  when  the  chord  length  is  finite  but 
acoustically  non-compact.  Howe’s  (2001a)  normalized  Green  function,  Gj,  shall  be  used  to 
model  this  effect.  This  yields  an  equation  for  the  dipole  sound  from  the  thick  but  acoustically 
non-compact  airfoil  in  terms  of  the  mean  velocity  and  the  up-wash  velocity  spectrum  and 
spanwise  integral  scale: 


The  normalized  Green  function  is  given  by 


Gj.(x,y,co)  =  l+ 


G,j;(x,y,a))  +  GrpXx,y,Q)) 

^lUplane  (-^>^>^0) 


(3.9) 


where  G^^,  Gje,  G^^piane  represent  respectively,  the  contributions  from  the  leading  edge,  the 
trailing  edge  and  the  governing  Green  function  for  the  “half-  slab”  diffracting  surface  used  in  Eq. 
(3.7).  These  are  given  by  Howe  (2(X)1)  with  [Fig.  (3.1)]  &=\i/=7Zf2: 
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(3.11) 


GiE{x,y,(0) 


iK^'^  U I  [1 + /(2mV)]  ^ 


=y(^)  +igix)  are  the  Fresnel  integral  auxiliary  functions  whose  approximations  are  given  by: 


1  +  0.926X  .  f  ■)_  1 

^  2  +  1.792x  +  3.104x'’^^^^"2+4.142x+3.492xV6.67a:'  ' 


(3.12) 


The  predicted  dipole  sound  from  Eq.  (3.8)  shall  be  compared  with  the  measured  sound 
reported  by  Paterson  and  Amiet  (1976)  in  Fig.  1.1.  Accordingly,  the  reported  up-wash  spectrum 
will  use  the  von  Karman  formulation  for  isotropic  turbulence  with  an  integral  scale,  ^  0.03  m, 
and  a  turbulence  intensity,  IIP  =  0.045, 


=  - IT 


(3.13) 


The  spanwise  integral  scale,  will  also  use  the  von  Karman  formulation. 


/3(®)  =  r 


_  sf  r(i/3)? 

3Lr(5/6)J  ' 


(4A,fc,/3)^ 

[3+ 8(4  A,*,  /3)']Jl+(4A.jt,/3)' 


(3.14) 


yielding  the  same  result  as  the  Corcos  model  or  the  Gaussian  model  for  frequencies  above 
300  Hz.  r  is  the  Ganuna  function  with  r(l/3)=2.68,  and  r(5/6)=1.13. 

The  airfoil  chord  length  is  c=  0.23  m,  the  span  is  L=0.53  m.  The  thickness  (h)  to  chord 
(c)  ratio  of  the  NACA  foil  is  hlc=0.\2.  The  dipole  sound  in  Fig.  1.1  was  measured  at  a  distance 
j^=2.25  m  from  the  mid-chord  of  the  foil.  The  Mach  number  correction  in  Eq.  (3.2)  becomes, 
Mo/e=0.05  M,  and  e=\if=nl2.  In  Fig.  3.3  for  Ink  =0  and  in  Fig.  3.4  for  hk  =0.12,  the  predicted 
dipole  sound  spectra  are  plotted  using  Eq.  (3.8)  at  velocities  of  40m/s,  60  m/s,  90  m/s  and  120 
m/s.  These  predictions  are  compared  at  each  velocity  to  the  measured  dipole  sound  from 
Paterson  and  Amiet  (1976, 1977). 


22 


120  m/s;  h/c=( 


o  o  o 

IL  If  I. 


^  ^  ^  B  S  B  3 

S  IS  is  ^  ^  T3  OJ 

jo  jo  jo  M  M  JO  > 

I  1  1  1  i  6  O 

g  s  §  ?  S  g  2 

IllUttl 


ZH/Bd  oJOjLU  oz  SJ  9P  ‘Ids 


23 


Figure  3.3  A  comparison  of  the  NACA  0012  (h/c=0A2)  dipole  sound  measurements  of  Paterson  and  Amiet  (1976, 1977) 
(points),  and  the  thickness  theory  [Eq.  (3.8)]  with  h/c=0.0  (line)  for  free  stream  velocities  of  40  m/s,  60  m/s,  90 
m/s,  and  120  m/s. 


h/c=0.t2 


Figure  3.4  A  comparison  of  the  NACA  0012  (A/c=0.12)  dipole  sound  measxirements  of  Paterson  and  Amiet  (1976, 1977) 
(points),  and  the  thickness  theory  [Eq.  (3.8)]  with  /i/c=0.12  (line)  for  free  stream  velocities  of  40  m/s,  60  m/s,  90 
m/s,  and  120  m/s. 


4.  Dipole  Sound  From  a  Thick,  Symmetric  Foil  in  a  Turbulent  Flow  With  Mean  Shear. 

In  Section  1,  the  turbulence  quadrupole  source  was  specified  to  be  without  mean  shear. 
These  types  of  turbulent  flows  are  typically  generated  by  turbulent  wake  flows  and  jet  induced 
shear  layers  that  have  had  sufficient  time  to  diffuse  or  by  the  outer  portion  of  a  turbulent 
boundary  layer  (TBL)  flow.  Conversely,  wake  and  jet  flows  have  regions  of  high  mean  shear 
before  sufficient  diffusion  occurs,  as  do  turbulent  flows  near  their  bounding  surfaces.  The 
portion  of  the  acoustic  source  that  contains  the  mean  shear-turbulent  (MS-T)  source  has  been 
used  to  describe  the  quadrapole  pressure  radiation  from  turbulent  shear  flows  that  occur  in  TBL 
flows  [Kraichnan  (1956),  Chase  (1980),  Blake  (1986)]  and  jets  [Moon  (1975)].  In  this  section, 
the  dipole  sound  due  to  a  MS-T  source  incident  to  the  leading  edge  of  a  thick  foil  will  be 
modeled. 

Referring  to  Fig.  4.1(a),  a  foil  of  thickness,  h,  and  chord  length,  c,  encounters  a  turbulent 
shear  layer,  dUi(y2)/dy2=const.,  where  the  shear  layer  is  thicker  than  the  foil  thickness.  The  force 
induced  by  the  diffraction  of  the  turbulent  shear  layer  by  the  leading  edge  is  formulated  as  a 
Kirchhoff  integral  for  the  scattering  of  the  incident  surface  pressure  by  a  “semi-infinite  slab” 
with  thickness,  h,  and  rounded  leading  edge  [Fig.  4.1(b)].  The  modeling  approach  is  identical  to 
that  described  in  Sections  2  and  3.  In  this  section,  the  wall  pressure  spectram  imposed  on  a  rigid 
plane  will  be  calculated  from  the  quadrapole  radiation  of  the  aero-acoustic  source  that  has  been 
specialized  to  be  linear  in  fluctuating  terms  and  to  contain  mean  shear.  This  pressure  spectram 
can  be  used  with  the  results  of  Section  2  to  calculate  the  leading  edge  noise  for  a  thick  foil 
cutting  through  a  shear  layer  that  will  be  specialized  to  be  spatially  homogeneous. 

At  low  Mach  numbers,  the  acoustic  pressure,  p,  generated  by  the  mean  shear-turbulence 
(MS-T)  source  is  described  by  the  following  acoustic  analogy: 

(4.1) 

cl  dr  dy2  dy^ 
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Figure  4.1  (b)  Representation  of  foil  as  a  semi-infinite  slab,  with  thickness,  h,  and  rounded  leading 
edge  for  leading  edge  force  estimate  due  to  diffraction  of  incident  surface  pressure,  p. 


The  pressure,  p,  and  the  fluctuating  up-wash  velocity,  w'a ,  are  in  the  time  domain  and  are  denoted 
by  bold  face  type.  The  planar  wave  number  spectrum  of  the  wall  pressure  that  is  radiated  by  the 
acoustic  source  in  Eq.  (4. 1)  becomes 


-ikx+ikar 


where  the  pressure,  p,  and  up-wash  velocity,  u^,  have  been  Fourier  transformed  to  the  frequency 
domain  and  are  no  longer  in  bold  face  type. 

The  planar  Fourier  transform  of  the  Green  Function  is  given  by 


1  ff 

—  ff  3  =  i- - 

27t^i  r  y(k) 


where 


'y(k)  =  Jk„-k^  -^3 


and  —  coICq,  >  yi>3—  (•^i>'^3) » <^Xi,'j—dxidxy  Substituting  Eq.  (4.3)  into  the 

Eq.  (4.2), 


Pik,a>)=-^  \\\dy^dy^^^ 


dU,  du2  ,  ..e'- 

«U  >1.30/1 - L - fv  11-  V-  — 


and  rearranging  the  terms  of  the  integrand  gives  the  following  expression  for  the  wall  pressure 
wave  number  spectrum: 


^2 
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The  Fourier  transform  of  the  fluctuating  up-wash  velocity  with  respect  to  the  in-plane 
coordinates  is  denoted  by  an  over  bar 


TW 

-«  [27t) 


and  Fourier  transform  of  the  derivative  of  the  up-wash  velocity  is: 


(iTT) 


Substituting  Eq.  (4.8)  into  Eq.  (4.6)  determines  the  wall  pressure  wavenumber  spectrum 


f  dU  — 

p(k,a))=  J  dy22p-^'——ik,U2(k„y2,k3,a>) 

yj=0  ^2  JW 


and  the  complex  conjugate  of  the  wall  pressure  spectrum  becomes: 


»  ■  ■  f  •  dU\  -r  ,  . 

pik,(0)=  J  dyilp-^i-i)  iriK)^i{kx,y 2,1c i,(6)  . 

,j=o  <^2  y{k) 


(4.10) 


The  wall  pressure  power  spectrum  is  then  determined  from  the  ensemble  average  of  the  pressure 
and  its  complex  conjugate: 


(k,ct))5(k  -  k  )S((0  -(»)=<  p(k,  (o)  p*  (Jc,(6)  > 
yielding  the  wall  pressure  wave  number  frequency  spectrum : 


(4.11) 


=  4p' J  rfy  J  dy2^?-^  (^,3.y2.3'2.®)l 


,()'2r(t)-o'2r  (*) 


^2  ^2 


.  (4.12) 


The  following  separable  form  to  the  up-wash  wave  number  spectrum, 


^u,uM^y2>y2>ki,(0)  =  <D(ti))0i(*,)i?22(3'2>3'2)03(^3) 


(4.13) 


will  facilitate  the  integration  of  Eq.  (4. 12).  The  wave  number  spectra  are  normalized 


ldkj(t>jikj)  =  i 


(4.14) 


and  the  streamwise  (1)  and  spanwise  (3)  wave  number  spectra  will  use  Corcos’  (1963)  and  a 
Gaussian  form  respectively: 


(4.15) 


The  shear  layer  will  be  considered  to  be  much  thicker  than  the  airfoil.  In  this  context, 
this  model  will  assume  that  the  surface  normal  correlation  length,  12=1.417^ /co ,  is  much  smaller 
than  the  thickness  of  the  shear  layer.  The  surface  normal  correlation  function  will  accordingly  be 
taken  to  be  spatially  homogeneous  with  correlation  length  4: 


A2(y2.)'2>®)  =  ^(3’2->’2)4(®)  • 


(4.16) 


It  is  specified  that  the  mean  shear  does  not  change  significantly  over  the  correlation 
length,  k.  Upon  insertion  of  Eqs.  (4.13),  (4.15)  and  (4.16)  into  Eq.  (4.12)  and  integrating  with 
respect  to  y2>  th®  wall  pressure  wavenumber  spectrum  becomes: 


(fe,_3,m)  =  4p^  J  ^>2^1^*^^—  (<»)A(^)<^3(^3)4(®) 


Jy2Y{k)-iy2r*(f^) 


(4.17) 
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The  wall  pressure  frequency  spectrum  is  determined  by  integrating  over  Jk,  and  ky 
Integration  over  the  streamwise  wave  number,  is  facilitated  by  assuming  that  the  correlation 
function  in  the  streamwise  direction  is  peaked  at  the  convection  ridge,  ik,  =k^.  The  streamwise 
wavenumber  spectrum  can  be  approximated  by  a  delta  function 


(4.18) 


yielding 


0  y  ^2  J 


(4.19) 


where 


TiK)  =  V*o  -  -  *3  ~  +*3  . 


(4.20) 


Substitution  of  Eq.  (4.20)  into  (4.19)  gives  the  spanwise  wall  pressure  wave  number  spectrum 


^,^{k„G))  =  4p^  J  O  {co^,{kM(o) 


-y2l2^k/+kl  ] 

V  +  ^3 


(4.21) 


The  spanwise  wave  number  spectrum  is  peaked  at  k^^O.  This  allows  the  spanwise  wave  number 
correlation  coefficient  to  be  approximated  by  a  delta  function: 


03(^3)  “^*3)- 


(4.22) 


Insertion  of  Eq.  (4.22)  into  (4.21)  yields  Eq.  (4.23): 


+»  /  \Z  +0O 

0  y  ^2  J  ^  k^  +k^ 


(4.23) 
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Integration  over  the  spanwise  wave  number  yields  the  wall  pressure  frequency  spectrum, 


=2p^  J 

0  V  ^2  y 


(4.24) 


and  integrating  with  respect  to  yields  the  point  pressure  spectrum  for  the  mean  shear  source: 


/  Y 

\  <^2  J 


(4.25) 


The  ratio  of  the  point  pressure  spectrum  from  the  mean  shear-turbulence  (MS-T)  source 
[Eq.  (4.25)]  to  the  source  without  mean  shear  (T-T)  [Eq.  (2.39)]  is: 


_  \du,  1 


^p/(®)  4.6m 


(4.26) 


This  defines  a  cut-off  frequency  below  which  the  dipole  sound  from  the  MS-T  source  should 
dominate  the  dipole  sound  from  the  source  without  mean  shear  (T-T). 

The  leading  edge  dipole  sound  is  determined  by  multiplying  the  ratio  of  the  MS-T  wall 
pressure  spectrum  to  the  T-T  spectrum  [Eq.  (4.26)]  by  the  dipole  sound  due  to  the  T-T  source 
[Eq.  (3.8)]: 


^radMs-T.  s  „  I G^. (fc,c sin(vr))  P  mLsin^(0/2)sinvr  tt  . 

^  il+ n\\x\^  K 


(  Y 

.  6.4p't/>„^„^(m)/3(m)/(;r)  . 

\dy^  4.6C0 J 


(4.27) 


The  dipole  sound  measurements  made  by  Olsen  and  Wagner  (1982)  with  airfoils  in  the 
shear  layer  of  a  10  cm  diameter  round  jet  will  be  predicted  using  estimates  of  the  incident 
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sources  with  and  without  mean  shear  from  equations  [Eq.  (4.27)]  and  [Eq.  (3.8)]  respectively. 
The  chord  length,  c,  of  the  foil  was  2.54  cm.  It  will  be  assumed  that  the  wetted  span  of  the 
airfoil,  L,  was  equal  to  the  jet  diameter.  The  airfoil  thickness  to  chord  ratio  ranged  from 
0.032  <h/c<  0.375.  The  flow  speed  incident  to  the  leading  edge  of  the  foil  was  94  m/s  and  the 
integral  scale  was  /y=1.6  cm.  The  fluctuating  velocity  was  25%  of  the  mean  velocity  incident  to 
the  leading  edge  of  the  foil.  The  mean  shear,  1723/sec  was  estimated  from  the  mean 

velocity  measurements  of  Wygnanski  and  Fiedler  (1965)  for  a  round  self-similar  jet.  The  sound 
pressure  levels  in  one-third  octave  frequency  bandwidths  were  measured  at  roughly  y2  =  4.57  m 
with  &=\ifi=jt/2  from  Fig.  3,1.  The  leading  edge  of  the  foil  was  four  jet  diameters  downstream  of 
the  lip  of  the  jet. 

The  predicted  narrow  band  sound  pressure  levels  in  dB  re  20  micro  Pa/Hz  for  both 
sources  are  plotted  in  Fig,  4.2  for  a  flow  speed  of  94  m/s.  The  predicted  dipole  sound  spectrum 
from  the  mean  shear  source  dominates  below  lOOHz.  The  dipole  sound  without  the  mean  shear 
source  is  predicted  for  foil  thickness-to-chord  ratios  of  h/c=  0.032  and  0.375.  Above  1  kHz,  the 
thickness  effect  starts  to  attenuate  the  predicted  dipole  sound. 

Figure  4.3  shows  the  measured  dipole  sound  from  Olsen  and  Wagner  (1982,  Fig.  2a)  in 
one  -third  octave  bandwidths.  The  thickness  theory  accounts  for  the  relative  differences  in  the 
dipole  sound  for  the  foils  as  a  function  of  thickness.  The  dipole  sound  is  not  measured  below 
200  Hz  where  the  predicted  mean  shear  dipole  sound  should  dominate.  This  may  explain  why 
the  theory  used  by  Olsen  and  Wagner  (1982)  was  able  to  predict  the  leading  edge  noise  for  the 
thin  airfoil  without  taking  into  account  the  mean  shear  source.  The  predictions  made  in  Fig.  4.2 
violate  the  “rapid  distortion”  approximation  since  the  turbulence  levels  are  above  10%.  The 
mean  shear  levels  were  estimated  assuming  self-similarity  of  the  mean  velocity  profiles  even  for 
distances  downstream  of  the  jet  exit  as  close  as  four  jet  diameters. 
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Predicted  leading  edge  dipole  sound  pressiire,  dB  re  20  micro  Pa/Hz  at  94  m/s  with  mean  shear  source 
/i/c=0.032  (blue)  and  without  mean  shear  (no  m/s)  for  /i/c-0.032  (red),  0.375  (green). 


AIRFOIL-AIONE  NOISE.  SPl^,  dB  (ref.  2k10’5  HImh 


5.  Conclusions 


The  dipole  sound  due  to  the  diffraction  of  homogeneous  turbulence  by  the  leading  edge 
of  a  thick  foil  was  calculated  using  the  Kirchhoff  surface  integral  with  a  rigid  surface  Green 
function.  The  Green  function  accounted  for  the  foil  thickness  and  the  acoustic  back  scattering 
associated  with  the  acoustically  non-compact  chord  length.  The  surface  stagnation  enthalpy  was 
calculated  with  a  volume  integral  incorporating  the  rigid  plane  Green  function  and  an  equivalent 
form  of  the  lifting  dipole  component  of  Howe’s  (1975)  acoustic  source  that  was  derived  in  terms 
of  the  mean  free  stream  velocity  and  the  fluctuating  up-wash  velocity  of  Sears’  (1941)  analysis. 
The  derivation  of  the  equivalent  form  of  the  source  required  that  the  velocity  field  be  solenoidal. 

An  estimate  of  the  dipole  radiation  was  made  using  the  foil  geometry  and  measured  grid¬ 
generated  turbulence  statistics  from  Paterson  and  Amiet  (1976, 1977).  The  predictions  show 
good  agreement  with  the  measured  dipole  sound.  This  establishes  that  the  square  of  the 
magnitude  of  Sears’  function,  for  a  foil  of  thickness,  h,  can  be  expressed  as 

l5(k„/i)P=IS(k,)Pexp[-ki/i/2],  (5.1) 

where  S(ki)  is  Sears’  function  with  h=0. 

The  modeling  approach  requires  that  the  incident  sources  obey  the  “rapid  distortion” 
approximation.  This  requirement  may  be  violated  when  the  turbulence  fluctuations  are  too  large 
relative  to  the  free  stream  geometry.  When  this  occurs,  the  sources  may  interact  with  each  other 
or  their  images.  The  RMS  turbulence  fluctuations  of  the  validation  data  of  Paterson  and  Amiet 
(1976, 1977)  were  approximately  4.5%  of  the  free  stream.  This  is  well  within  the  limits  (less 
than  10%)  established  by  Grace  (2001)  to  ensure  that  the  sources  do  not  deviate  from  the 
potential  flow.  Leading  edge  flow  separation  from  a  sharp  edge  may  also  cause  the  trajectories 
of  the  incident  sources  to  deviate  sufficiently  from  that  modeled  by  the  Green  function.  This 
modeling  approach  also  required  the  use  of  the  frozen  flow  approximation  that  is  strictly 
applicable  to  constant  velocity  flows  in  the  vicinity  of  the  diffraction  zones.  The  success  of  the 
dipole  sound  predictions  suggests  that  the  frozen  flow  approximation  may  be  applied  to  steady 
flows  whose  potential  flow  streamlines  exhibit  curvature  and  consequential  acceleration  of  the 
sources  around  an  edge.  Extensions  of  the  modeling  approach  can  be  made  to  predict  the  affect 
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of  foil  camber  on  the  dipole  sound  by  incorporating  this  effect  into  the  velocity  potential  of  the 
Green  function  via  a  conformal  mapping,  as  demonstrated  by  Martinez  and  Rudzynsky  (1997) 
for  the  incompressible  problem. 

The  dipole  sound  from  a  foil  cutting  through  a  mean  shear  layer  was  obtained  using  the 
acoustic  source  that  was  tailored  for  mean  shear  sources.  The  incident  surface  pressure  spectrum 
due  to  this  source  that  is  imposed  on  a  rigid  plane  was  calculated.  The  ratio  of  the  dipole  sound 
of  the  mean  shear  source  to  the  source  without  mean  shear  was  determined  to  be  proportional  to 
the  square  of  the  ratio  of  the  mean  shear  to  the  frequency  of  the  sources.  This  establishes  a 
critical  frequency  below  which  the  dipole  sound  from  the  mean  shear  should  dominate  and  above 
which  the  dipole  sound  from  the  source  without  mean  shear  should  dominate.  This  result  has 
been  reported  for  the  frequency  dependence  of  the  mean  shear  and  non-mean  shear  quadrupole 
sources  associated  with  jet  radiation  (e.g.,  Moon  (1975)). 

Estimates  of  the  measured  leading  edge  dipole  sound  by  Olsen  and  Wagner  (1982)  were 
made  with  and  without  the  mean  shear  incident  sources.  The  frequency  below  which  the  mean 
shear  source  was  estimated  to  dominate  was  lOOHz.  The  measured  dipole  sound  was  made  for 
frequencies  above  200  Hz.  This  may  explain  why  the  predicted  dipole  sound  made  by  Olsen  and 
Wagner  (1982)  agreed  with  the  measurements  even  though  they  did  not  take  into  account  the 
mean  shear  source. 
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